Rare earth spin ensemble magnetically coupled to a superconducting resonator 
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Interfacing superconducting quantum processors, working in the GHz frequency range, with opti- 
cal quantum networks and atomic qubits is a challenging task for the implementation of distributed 
quantum information processing as well as for quantum communication. Using spin ensembles of 
rare earth ions provides an excellent opportunity to bridge microwave and optical domains at the 
quantum level. In this letter, we demonstrate magnetic coupling of Er 3+ spins doped in a Y2SiOs 
crystal to a high-Q coplanar superconducting resonator. 



Quantum communication is a rapidly developing field 
of science and technology which allows the transmission 
of information in an intrinsically secure way [1 j. As well 
as its classical counterpart, a quantum communication 
network can combine various types of systems which 
transmit, receive, and process information using quan- 
tum algorithms [2]. For example, the nodes of such net- 
work can be implemented by superconducting (SC) quan- 
tum circuits operated in the GHz frequency range [3], 
whereas fiber optics operated at near-infrared can be used 
to link them over long distances. For reversible transfer 
of quantum states between systems operating at GHz 
and optical frequency ranges one must use a hybrid sys- 
tem [4]. Spin ensembles coupled to microwave resonator 
represents one of the possible implementations of such a 
system [SJ [6] . The collective coupling strength of a spin 
ensemble is increased with respect to a single spin by the 
square root of the number of spins. Transparent crystals 
doped with paramagnetic ions often possess long coher- 
ence times [7]-[9], and the collective coupling has been 
recently demonstrated with NV-centers in diamond [TOl - 
[12] and (Cr 3+ ) ions in ruby [Up. 

In this letter, we report on the measurement of a spin 
ensemble of Erbium ions in a crystal, magnetically cou- 
pled to a high-Q coplanar SC resonator. The Er 3+ ions 
are distinct from other spin ensembles due to their opti- 
cal transition at telecom C-band, i.e. inside the so-called 
"Erbium window" at 1.54 jiva wavelength, and long mea- 
sured optical coherence time. 

The energy level diagram of Erbium ions embedded 
inside a crystal is shown on FigjTJa). The electronic 
configuration of a free Er 3+ ion is 4/ 11 with a 4 I term. 
The spin-orbit coupling splits it into several fine struc- 
ture levels. An optical transition at telecom wavelength 
occurs between the ground state 2s+1 Lj = 4 I 15 / 2 and the 
first excited state 4 Ii3/2, where S, L and J are spin, or- 
bital and total magnetic momenta of the ion. The weak 
crystal field splits the ground state into eight (J+l/2) 
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FIG. 1: (Color online) (a) The structure of the energy levels 
of Er 3+ ions, (b) A picture of the rare-earth ion chip. The 
Er:Y2Si05 crystal of 1 x 1.5 x 3 mm 3 size is placed on the 
A/2 resonator. The width of the central transmission line of 
the resonator is 20 /xm; gaps are 9 /im. The magnetic field is 
applied along the chip surface, (c) The sketch of the crystal 
orientation marked with its optical extinction axis b and Di 
with respect to directions of the bias magnetic field B and 
microwave field B(t). 



Kramers doublets [15] . At cryogenic temperature, only 
the lowest doublet Zi is populated, therefore the sys- 
tem can be described as an effective electronic spin 
with S — 1/2. However, Erbium has five even isotopes 
162 Er, 164 Er, 166 Er, 168 Er and 170 Er, and one odd isotope 
167 Er (natural abundance 22.9%) with a nuclear spin 
1=7/2. Therefore, the electronic states of 167 Er with ef- 
fective spin projection ms = ±1/2 are additionally split 
into 8 hyperfine levels [T4] , 

The magnetic properties of Erbium ions are associ- 
ated with an unquenched total orbital moment J in a 
crystal field which results in the appearance of a large 
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magnetic moment of nearly 7/i#, where \±b is the Bohr 
magneton, at particular orientations of the applied mag- 
netic field [15]. The large spin tuning rate ~ 200 GHz/T 
makes Erbium doped crystals favorable for its integra- 
tion with SC qubits that can be operated only at rel- 
atively low magnetic field. Such an integrated spin-SC 
device can be very attractive for its applications particu- 
larly in quantum repeaters, where one can store quantum 
information and perform local operations [16 j. A quan- 
tum state of the optical field can be mapped into spin- 
waves encoded in Zeeman or hyperfine levels of Erbium 
ions [U] and transferred later to the quantum state of the 
microwave field [4] [6]. The strong coupling between SC- 
qubits and microwave resonator allows fast quantum gate 
operations on the state of microwave field at nanosecond 
timescales [THJ [19] . This sequence can also be launched 
in reverse order, thus establishing a coherent quantum 
transfer between GHz and optical frequency ranges. 

In this experiment, we use a single Y 2 Si05 crystal 
doped with 0.02% of Er 3 + (Er:YSO), supplied by Sci- 
entific Materials Inc. The crystal has dimensions of 
1 x 1.5 x 3 mm 3 and it is glued on top of the sili- 
con chip with a A/2-coplanar Niobium SC-resonator, see 
Figjljb). The resonance frequency of the rare-earth ion 
chip is uo = 27rx8.9 GHz and its quality factor is Q w 
32000 at magnetic field of Erbium transitions. The crys- 
tal orientation is shown in Fig[TJc) with its optical ex- 
tinction axes b and Di and is specified by the angles 
and <p between these axes and direction of the applied 
magnetic field [20] . The b-axis of the crystal is directed 
along its 1.5 mm side and is perpendicular to the bias 
magnetic field B (9 = 90°) applied parallel to the chip 
surface. The angle between the magnetic field and Di 
is (j) = —60°. The particular orientation of our experi- 
ment is chosen to maximize the g-factor for the crystallo- 
graphic site 1 [14] [20] , which falls into our field scanning 
range between and 70 mT and is relevant for our experi- 
ment. The gap between the bottom of the crystal and the 
chip surface is controlled by observing Newton's interfer- 
ence fringes and with proper placement 1-2 fringes are 
visible, yielding a gap of about 0.5 fim. The experiment 
was performed in a 3 He cryostat, at a base temperature 
of 280 mK. 

In quantum optics, the resonator- spin ensemble in- 
teraction is usually described by the Tavis-Cummings 
Hamiltonian [21 . Provided that number of photons in 
the resonator is much smaller than number of spins TV, 
a spin ensemble behaves as a harmonic oscillator cou- 
pled to a cavity. That results in an avoided level crossing 
when the spins are tuned into resonance with a cavity, 
which is also in agreement with a phenomenological treat- 
ment of the observed effect. The rotating component of 
the magnetization of the spins produces and oscillatory 
magnetic field, which perturbs the inductance Lq of the 
resonator. The new inductance V = Lo(l + x£), where 
X = x' ~ ix" is the dynamic magnetic susceptibility of 



the ensemble with x' and x" being its dispersive and ab- 
sorptive parts respectively, and £ is a geometric factor 
describing the spin distribution across the mode and its 
coupling to the oscillating field. Thus, the resonance fre- 
quency of the rare-earth chip is uj'q = ojo/ \/l + x£- P ro_ 
vided that x£ we obtain the following equations for 
the resonator frequency u f and its decay rate ^q: 

/ , V<2 (^0 -U a ) 

v 2 r* 

K ° = Ko+ K- Ws )2 2 +rf ' (2) 

where v is coupling strength, uj s is the Larmor frequency 
of spin ensemble, and is its total linewidth [13]. The 
coupling strength v can be expressed via the static sus- 
ceptibility of the spin ensemble xo and reads as v = 
ujo \/xo£/2. The substitution of the actual value of xo 
yields the final expression for the collective coupling 
strength v = gfiB v^¥^/^? where g is an effective 
g-factor due to the magnetic anisotropy of the crystal 
and n is the spin concentration. For the concentration of 
Erbium electronic spins of ns ~10 18 cm -3 , the effective 
g-factor in the plane perpendicular to the magnetic field 
is g ~ 7 and filling factor £ ~ 0.25 we expect the coupling 
strength v/2ir - 60 MHz. 

The microwave transmission spectroscopy of the rare- 
earth ion chip as a function of the bias magnetic field 
is presented in Fig[2] By using a vector network an- 
alyzer we measured the S21 parameter, which contains 
both the amplitude and the phase of the signal. The 
probing power at the input of the resonator is about 1 
fW corresponding to an excitation level of ~ 100 mi- 
crowave photons in the SC-resonator. The amplitude of 
the transmitted signal IS21 1 as a function of magnetic field 
and probing frequency is presented in Figj2|a). The spec- 
trum consists of the large avoided level crossing at a bias 
field of B = 55 mT and associated with magnetic dipole 
transition between electronic states with ms = ±1/2. 
A regular pattern of six dark interruptions is associated 
with the odd Erbium isotope and is due to allowed hy- 
perfine transitions between states with equal nuclear spin 
projection mi = 7/2, 5/2,..., -1/2, -3/2. The states with 
mi = -5/2 and -7/2 do not enter into our scanning range. 

To understand the measured spectrum we numerically 
diagonalize the spin Hamiltonian for paramagnetic ion in 
the crystal [T3lfT4] : 

H = fi B B gS + IAS + IQI- /i n g n B • I, (3) 

where g is the g-factor tensor, A is the hyperfine tensor, 
Q is the nuclear quadrupole tensor, \i n is nuclear Bohr 
magneton and g n is nuclear g-factor. The first term in 
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FIG. 2: (Color online) Transmission spectroscopy of the rare- 
earth ion chip, (a) The amplitude of the transmitted signal 
| S21 1 • The energy spectrum of different magnetic transitions is 
shown with dashed lines. The spin tuning rate js/^ty = 162.9 
GHz/T and values for the hyperfine splitting are taken from 
the fit of the frequency shift (see the text for details), (b) The 
calculated energy spectrum of Er:YSO in the magnetic field. 
Allowed and forbidden transition are shown with arrows. 

the Hamiltonian presents an electronic Zeeman splitting, 
the second one describes hyperfine interactions, the third 
one is the quadrupole term and the last one is the Zee- 
man splitting due to the nuclear spin. The values for 
the tensors has been taken from a previous ESR study 
of Er:YSO crystal [14]. The resulting eigenspectrum is 
presented in Figj2^c). The position of each energy level is 
drawn as a function of the applied magnetic field, and the 
experimentally observed magnetic transitions are shown 
by arrows. For the even Erbium isotopes only the first 
term in the Hamiltonian (|3| survives and that results in 
the strong magnetic transition between electronic spin 
states ms = ±1/2 shown with the red arrow and marked 
by the letter "S". The odd 167 Er ion has 8 allowed hy- 
perfine transitions when Ami = 0, which are shown by 
blue arrows and marked with m/ numbers. 

To extract the coupling strengths and linewidths for 
different transitions, each spectral line |S2i| is fitted to a 
Lorentzian at every value of the magnetic field. The data 
corresponding to the shift of the resonator frequency uj' 
for transitions "7/2" and "S" is shown in Figj3ja),(b). 
The dispersive behavior of the frequency shift in the 



vicinity of Erbium spin transitions is well fit with Eg.Q. 
For the transition "7/2", the coupling strength extracted 
from that fit is Vj^/^it = 2.1 ± 0.3 MHz, which ex- 
ceeds the decay rate of the resonator k/2tt — 0.3 MHz. 
However the linewidth of spin ensemble exceeds the cou- 
pling strength and the same fit yields Ty 2 /2ir = 70 ± 1 
MHz. The frequency shift of the resonator due to the 
coupling to the electronic spin ensemble "S" ( FigJSJb)) 
has been studied for two excitation levels of 10 2 (light 
gray dots) and 10 5 (dark gray dots) microwave photons. 
The dashed line described by Eq.Q fits well the low ex- 
citation spectrum yielding vs/2tt = 11.6 ± 0.2 MHz and 
T^/27T = 145 ± 5 MHz. At high excitation level, the pre- 
sented data corresponds to a Dysonian line shape associ- 
ated with spin diffusion in bulk metals [22]. We believe 
that the observed mixture of the dispersive and absorp- 
tive line shapes might also be associated with spin diffu- 
sion in and out of the part interacting with the resonator 
mode. 
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FIG. 3: (Color online) Frequency shift of the resonator 
uj'q /2tt— 8.9 GHz and its quality factor Q versus the applied 
magnetic field, (a) For the transition "7/2". Gray dots are 
experimental data. The dashed line is fit to Eq.|l]). (b) For 
the transition " S" . Light gray and dark gray dots are experi- 
mental data taken at excitation level of 10 2 and 10 5 photons 
respectively. The dashed line is a fit of the low excitation 
data to Eq.|l]). (c) Gray dots are experimental data points. 
The solid line is fit of experimental data. Clearly observed 
dips corresponds to different magnetic transition. Inset: The 
temperature dependence of the collective coupling v s versus 
temperature. Circles present measured data. The solid line 
is a fit to the data (see the text for details). 

The measured spectrum in FigJ5|a) also reveals other 
interesting features: between hyperfine transition one can 
recognize an additional regular structure. To study this, 
we plot the quality factor of the resonator on the Figj3^c) 
as a function of magnetic field. The curve consists of a se- 
ries of regular absorption dips originating from the mag- 
netic coupling of electronic and hyperfine spin ensembles. 
It also contains an additional weak pattern appearing in 
between the hyperfine transitions. We interpret these 
weak absorption lines as corresponding to the forbidden 
quadrupole transitions between the hyperfine states sat- 
isfying Ami = ± 1 [l3j[T4]. These transitions are marked 
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by the letters "Ql", "Q2" and "Q3". 

The experimental data on the Fig (3] can also be fit with 
Q = cuo / k'o by using Eq. ([2| , where each magnetic transi- 
tions contributes independently to the cavity decay rate 
K,' . The data points between transitions "Ql" and "5/2" , 
and after transition "-1/2" are not fit well due to the pres- 
ence of additional magnetic transitions related the other 
crystallographic site. The fit of the quality factor behav- 
ior yields the following values: v 7 / 2 /2ir — 2.5 ± 0.3 MHz, 
r; /2 /27r = 65 ± 3 MHz and v s /2tt = 13.8 ± 0.1 MHz, 
T* s /2tt = 141 ± 3 MHz. The slight differences of fit- 
ted parameters by using u f or Q are probably due to an 
admixture of two magnetic classes into one anticrossing. 

The dimensionless parameter to identify the coupling 
regime is the cooperativity C = v 2 / kT* [23], which cor- 
responds to the number of coherent oscillations between 
coupled systems. The cooperativity parameter for the 
hyperfine transition "7/2" is measured to be C 7 / 2 = 0.4. 
High cooperativity coupling is reached for the electronic 
transition "S", where Cs = 5.2. In that regime mi- 
crowave photons inside the SC-cavity coherently interacts 
with a spin ensemble. However, to observe a normal- 
mode splitting the condition 2v > (ft, T*) has to be ful- 
filled. 

The coupling strength v$ of the electronic spin ensem- 
ble is found to vary as function of its temperature. The 
transmission spectrum of the rare-earth chip has been 
taken at temperatures of 0.3, 0.5, 1 and 1.5 Kelvin. The 
coupling strength vs is extracted from the change of the 
Q and plotted as a function of the temperature in the 
inset of Figj2^b). The temperature dependence is fit to 
v s = v s (0) [tanh(ftcj / 2 ^£ T )] 1/2 , see [13 , and plotted 
with solid line. The coupling strength of the ensemble 
at zero temperature, vs(0)/2tt = 17.3 ± 0.2 MHz, is the 
only fitted parameter. 

We also found that linewidths of different spin ensem- 
bles are not the same and grow with magnetic field. We 
explain this by a small misalignment of the crystal with 
respect to the magnetic field. When the angle deviates 
from 90°, the additional degeneracy due to C 2 symmetry 
of the crystal is lifted and each Erbium transition splits 
further into two different magnetic subclasses [I4j [T5J [20] . 
Using the "Easyspin" package [24] we simulated ESR 
spectra of the Er:YSO crystal and bounded the maximum 
misalignment of the crystal to be AO <0.4°. Assuming 
that the observed effect is dominated by the magnetic 
class splitting, an additional contribution to the spin 
linewidth is estimated to be about 30 MHz. Such large 
linewidth can not be completely attributed to the de- 
phasing of spins. We measured the pure dephasing time 
using a Hahn echo sequence in pulsed ESR spectrometer 
at the temperature of 7 K, and found that for the "7/2" 

(7 /2) 

transition T 2 ~ 540 ns and for the "S" transition 

(S) 

T 2 ~ 200 ns. Moreover, at the temperature of 0.3 K 
we would expect T 2 to be longer by at least an order of 



magnitude [HJ[9]. Therefore, the additional contribution 
to the spin linewidths may be associated with field inho- 
mogeneity or surface magnetism of the SC-resonator [25J. 
In conclusion, we have presented measurements on an 

Er 3+ 

spin ensemble in Y2Si05 crystal magnetically cou- 
pled to a SC coplanar resonator. This hybrid system is 
a promising candidate for an interface between SC quan- 
tum circuits and optical quantum networks. The mea- 
sured coupling strengths of different spin ensembles en- 
coded in Er:YSO crystal exceeds the decay rate of the 
cavity. The demonstrated on-chip ESR spectroscopy al- 
lows us to resolve all spin state transitions of the Er 3+ 
ions. The presented experiment realizes a first step to- 
wards implementation of a rare earth based quantum 
memory with telecom wavelength conversion that could 
constitute an essential building block of future quantum 
networks. 
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